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Abstract 



o 

^vq , The symbol error rate of the minimum distance detector for an arbitrary multi-dimensional constellation 



impaired by additive white Gaussian noise is characterized as the product of a completely monotone function 
with a non-negative power of the signal to noise ratio. This representation is also shown to apply to cases when 
the impairing noise is compound Gaussian. Using this general result, it is proved that the symbol error rate is 
completely monotone if the rank of its constellation matrix is either one or two. Further, a necessary and sufficient 
c/3 ' condition for the complete monotonicity of the symbol error rate of a constellation of any dimension is also 

obtained. Applications to stochastic ordering of wireless system performance are also discussed. 
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^ ' I. Introduction 

\o ' 

O 

An important performance metric in digital communications is the symbol error rate (SER). The 
convexity of the SER in the signal to noise ratio (SNR) plays a critical role in various optimization 
^ ' problems [0, B3). Convex SERs have a negative first derivative and a positive second derivative with 
respect to the SNR. If all the successive derivatives of the SER also alternate in sign (referred to 
as complete monotonicity), then it is possible to express the SER as a positive mixture of decaying 
exponentials, which has applications in SER analysis over fading channels, as described in @J. This 
serves as the motivation to explore the complete monotonicity (cm.) properties of the SER of arbitrary 
multi-dimensional constellations. An overview of the literature addressing related properties of the SER 
follows next. 
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One-dimensional and two-dimensional constellations have been adopted in many communication sys- 
tems in the literature, and investigations into the properties of the SER of these constellations have 
revealed the convexity of the SER with respect to the signal-to-noise ratio (SNR) under impairing additive 
white Gaussian noise (AWGN) [0, Q. Some special cases of two dimensional constellations such 
as M-ary phase shift keying (M-PSK) and M-ary quadrature amplitude modulation (M-QAM) have 
SERs which are known to be completely monotone functions of the SNR 01, J6]|, which is a stronger 
condition than convexity. On the other hand, constellations of dimensionality greater than two (which we 
refer to as "higher dimensional constellations" henceforth) have found practical applications in satellite 
communications 0, jH and more recently, in optical communications flU, OH- Investigations of the 
convexity properties of the SER of such constellations are relatively scarce in the literature. It is known 
that the second derivative of the SER of a constellation of dimensionality greater than two is non-negative 
at sufficiently high SNR [5J. Although this result is a general one, it does not provide a conclusive 
means of determining whether a given arbitrary constellation has a convex SER or not. For certain higher 
dimensional constellations, analytical expressions for the SER have been derived in the literature, which 
can be used to deduce the convexity of the SER. For example, the class of constellations of dimensionality 
2, 3, 4 and 5 described in [[Til can be verified to have convex SERs under AWGN, by differentiation of the 
analytical SER expressions given in [fTTTl . On the other hand, verifying if a SER is completely monotone 
is difficult, even if a closed form expression for the SER is available. It has been established recently, that 
if the rank of a constellation matrix is at most two, then it will have a completely monotone SER under 
AWGN with ML detection JT). However, there has been no investigation into the complete monotonicity 
properties of the SER for higher dimensional constellations, which can be used to simplify the expressions 
for average SERs over fading channels, and to establish useful comparisons of average SERs of a system 
under two different fading channels [@], using the tools of stochastic ordering, which is surveyed next. 

The theory of stochastic orders provides a comprehensive framework to compare two random variables 
(RVs) [fT2l . and can be applied to compare two fading distributions based on a range of performance metrics 
such as the ergodic capacity or the SER. One example is the convex order, which finds applications in 
comparing the variability of two RVs by comparing the expectation over all convex functions over the two 
different distributions that are being compared [4j. There are many other stochastic orders that capture 
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comparisons of RVs in terms of size and variability. These include the Laplace transform (LT) order [fl2l 
which compares the real valued Laplace transform of the probability density function (PDF) of two positive 
RVs. The instantaneous SNRs corresponding to many fading envelope distributions such as Nakagami and 
Ricean exhibit LT ordering with respect to their respective line of sight parameters [4J. This ordering can 
be systematically exploited through its connection with cm. functions, yielding generic comparisons of 
averages of a cm. function (such as SER) over two different positive random variables, even in cases 
where closed-form expressions for these averages are unavailable. 

In this work, it is shown that the SER of an arbitrary multi-dimensional constellation impaired by 
additive independent and identically distributed (i.i.d) Gaussian noise under maximum likelihood (ML) 
detection can be represented as a product of a completely monotone function of SNR, and a power of SNR, 
which depends only on the rank of the constellation matrix. This result also generalizes to SERs under 
compound Gaussian noise, which includes many non-Gaussian noise distributions such as Middleton 
class-A noise [TT3l and symmetric alpha-stable noise lfl4l . The SER of an arbitrary multi-dimensional 
constellation is shown to be completely monotone if the constellation matrix has a rank of one or two. 
Since complete monotonicity implies convexity, the SER is a convex function of the SNR, provided that 
the constellation matrix has a rank of one or two. For constellations matrices whose rank is greater than 
two, it is shown that the complete monotonicity of the SER depends on the constellation geometry and 
choice of prior probabilities. This work also describes a novel stochastic order for fading distributions, 
which can be used to order the average SERs of arbitrary multidimensional complex constellations over 
quasi-static fading channels, and generalizes the existing Laplace transform order on random variables. 

The rest of the paper is organized as follows: Section [I] surveys the relevant mathematical background, 
including a summary of stochastic ordering theory. Section [III] describes the system model. The result on 
the representation of the SER of an arbitrary multi-dimensional constellation is detailed in Section [IV] In 
Section [V] the applications such as (i) ordering the average SERs of constellations with cm. SERs over 
two different fading channels, and (ii) comparing the average SER of an arbitrary constellation using a 
newly proposed stochastic order, are discussed. Conclusions are presented in Section [VD 

Here are some remarks on the notations used in this paper. Vectors are denoted by boldface lower-case 
letters. Sets are denoted using upper-case script font. The set of real and natural numbers are denoted by ft 



and N respectively. Transpose of a vector is denoted by (•) T . The multivariate real (circularly symmetric 
complex) Gaussian distribution with mean vector a, and covariance matrix C is denoted by jV(a, C) 
(CM(a, C)). E [g'(X)] is used to denote the expectation of the function g(-) over the distribution of the 
random variable X. The identity matrix is denoted by I. P Cji (p) denotes the SER conditioned on the i th 
symbol of the constellation being transmitted. P e (p) denotes the SER averaged over the constellation. 
P e (p) denotes the average SER over a fading channel (where the averaging is over the constellation points 
and the fading channel statistics). ||-|| denotes the vector-2 norm for both real and complex vector spaces. 
The indicator function is defined as I[x G S] = 1, if x G § and 0, otherwise. 

II. Mathematical Preliminaries 

A. Complete Monotonicity 

A function g : (0, oo) — ^ is completely monotone (cm.), if and only if it has derivatives of all orders 
which satisfy 

( ~ ir d^ (x) " °' 0) 
for all n E NU{0}, where the derivative of order n = is defined as g(x) itself. The celebrated Bernstein's 
theorem [15J asserts that, g : (0, oo) — > H is cm. if and only if it can be written as a mixture of decaying 
exponentials: 

oo 

g(x) = J exp(— ux)p(u)du, (2) 

Q 

where p : [0, oo) — > [0, oo) is called as the representing function of g, in this paper. It is straightforward 
to verify that cm. functions are positive, decreasing and convex, and that positive linear combinations of 
cm. functions are also cm. 

A function g : (0, oo) — V % is said to be completely monotone of order a G K if and only if x a g(x) 
is cm.. If g is cm. of order a, then g is also cm. of order 0, where < < a. In lfT6l Theorem 1.3], 
it is shown that a necessary and sufficient condition for g to be cm. of order a is that g(x) must be 
represented in the form (0), where the integral converges for all x > 0. In addition, p must be a — 1 
times differentiable, with the k th derivative of p(u) equal to zero at u = for < k < a — 2, and 
d a_1 p(u) I 'du^ 1 nonnegative, right-continuous and non decreasing. 
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B. Stochastic Orders 

Let S denote a set of real valued functions g : Ji + — > Ji, and X and Y be non-negative RVs. The 
integral stochastic order with respect to a set of functions S is defined as IfTTl : 

X < s y & E [g{X)) < E [g(Y)} , Wg e 9 , (3) 

where S is known as a generator of the order < s . We write X ^ s y, if © does not hold. Two examples 
of integral stochastic orders relevant to this paper are now considered, by specifying the corresponding 
generator set of functions S- 

1 ) Convex Order: In this case, S is the set of all convex functions, and the order is denoted as X < cx Y. 
Applying © using g(x) = x and g(x) = —x (both of which are convex), we have E [X] = E [Y] whenever 
X and y are convex ordered. Therefore, convex ordering establishes that the RVs have the same mean, 
and X is "less variable" than Y. 

Instantaneous SERs of two-dimensional modulations over fading channels are known to be convex 
functions of the instantaneous SNR [5J. Therefore, establishing convex ordering of two RVs can enable 
comparisons of the average performance of the same system subject to fading modeled by two different 
RVs. 

2) Laplace Transform Order: This partial order compares random variables based on their Laplace 
transforms. In this case, S = {g{%) I g{x) = — exp (— px) , p > 0}, so that Xi < Lt X 2 if and only if 

E [exp(-pX 2 )] < E [exp(-pXi)] , V p > . (4) 

One useful property of LT ordered random variables, which can be obtained using © and © is that 
X\ <Lt X 2 is equivalent to 

E [g(X 2 )} < E \g{X 1 )] , (5) 

for all cm. functions g(-). In other words, the generator S can be replaced by the set containing the 
negative of all cm. functions without changing the stochastic order [[171 . In a wireless communications 
context, let p > represent the average SNR, and pX\ and pX 2 represent the instantaneous end-to- 
end SNRs of a system over two different fading channels with X\ <Lt X 2 , such as the case when 
y/X~i and y/X~2~ are Nakagami distributed with parameters mi and m 2 respectively, where mi < m 2 [ffl. 
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Let P e (px) be the instantaneous SER of a modulation scheme, and let P e (px) be cm.. For example, 
P e {px) = (1/2) exp(— px) for the case of differential PSK modulation, or P e {px) = Q (y/Zpx) for BPSK 
(where Q (y/x) : = /^(l/v 7 ^) exp(-t 2 /2)dt). Then using ©, it is seen that E [P e {pX 2 )\ < E [Pjj)X x )\, 
Vp > 0. This is because both choices of P c {px) are cm. functions for all p, and © is used with 
g(x) = P e (px). In other words, © can be used to compare the average SERs of the two channels at all 
SNR. 

C. Polyhedra and Polytopes 

A set CP C Ji N is a polyhedron if it is the intersection of finitely many closed half-spaces, i.e., CP := 
{x|Ax < b}, for some A G Ji MxN , h G Jl N , and the inequality is applied component- wise |[T8l . A 
bounded polyhedron is referred to as a polytope. In a digital communications context, it is known that the 
decision region of a multi-dimensional constellation impaired by white Gaussian noise is a polyhedron. 

A face of a polyhedron CP is the intersection of CP with a supporting hyperplane of CP. The dimension 
of a face is defined as the dimension of its affine hull. Faces of CP of dimension zero, dimension one, and 
dimension A — 1 are known as vertices, edges and facets of CP, respectively. 

Some examples of polyhedra relevant to this paper are described next. A polyhedral cone is a polyhedron, 
which is defined as cone(y) := {$^j=i A^y* | y* G V, Aj > 0}, where ^ is a non-empty set of points in 
Jl N . If the elements of V are linearly independent, then the polyhedral cone is called a simplicial cone. A 
well known result in combinatorial geometry literature is that any polyhedral cone admits a decomposition 
into simplicial cones |[P9l Lemma 1.40]. In this context, a decomposition of a polyhedron CP is defined 
as a collection of sets {CP 1; . . . , CP#}, such that U^CP,. = CP, and the intersection of any two sets in the 
decomposition is a common face of both, or the empty set. 

III. System Model 

In this paper, the transmission of an A-dimensional baseband signal s through AWGN is considered, 
which is described as follows: 

y = s + z , (6) 



where the transmitted symbol s G Jl N is chosen from § := {s!,...,s M }. The constellation matrix 
corresponding to S is defined as S := [s 1; . . . ,s M ], and the reduced dimension of S is defined as the 
rank of S, which is denoted by iV*. In ©, the noise is assumed to be z ~ A/"(0, (l/p)I), p > 0. When 
the signal energy is normalized as M _1 J2i=i ll s «l| 2 = 1> men Af -1 X]f=i ll s «l| 2 /^ [ zTz ] = P represents 
the average SNR per dimension. At the receiver, the ML detector under AWGN is assumed, where the 
detected symbol s is given by: 

A .. 1 1 2 

s = argmin | |y — s| | , (7) 

ses 

which is the ML detector for white Gaussian noise. Assuming that the origin of the coordinate system is 
shifted to Sj, the Voronoi region of Sj, denoted by %i is given by 

%i := {x G $. N \AiX < hi}, (8) 

where Aj G Ji FxN , where F < M, and the j th row of Aj is aj^ = (sj — Sj) T / ||s 3 - — Sj||, while the 
j th element of bj is bj^ — \\sj — Sj|| /2. It is assumed that © is a non-redundant description of DQ. The 
minimum distance d min of the constellation is defined through its square as <i min := min | |sj — Sj \ | 2 = 

Si ,Sj GS;Sj T^Sj 

4min tt.. Il20l. 

The probability of error P e i (p), given that s» is transmitted is given by 



p T 

-XX 

2 



dx , (9) 



where Si — S 2 : = {x G §i|x ^ S 2 } is the relative complement of S 2 in §> 1 . The probability of error 
averaged across all possible transmitted symbols is given by 

M 

P c (p) = ^Pr[s = s l ]P c ,,,(p) , (10) 
i=i 

where Pr[s = Sj] represents the a priori probability of transmitting Sj. 

For any constellation S, there exists an equivalent full rank constellation S*, which has the same SER 
as that of S. Such a definition is useful in developing a representation for the SER of a multidimensional 
constellation. 

Definition 1 (Reduced Constellation). Let S = USV T be a singular value decomposition of S, where 
U G R NxN , £ G Jl NxM and V G % MxM , and the diagonal matrix consisting of the first N* rows and 
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first N* columns of X contains the non-zero singular values of S. Then the N* x M matrix S* given by 
the first N* rows of S V T is defined as the reduced constellation of S. 

By definition, the number of rows of S* is no greater than that of S. In addition, S* can be shown 
to have the same SER as S. To see this, recall that the SER of the minimum distance detector depends 
only on the distance between the column vectors of S and the Frobenious norm of S. It then suffices 
to show that the columns of S* have the same distance properties and norm as that of S. To this end, 
observe that SV T = U T S has the same distance and norm properties as that of S, since it is an orthogonal 
transformation on S. Further, by construction, the last N — N* rows of SV T are zeros. Hence, the distance 
and norm properties of SV T and S* are identical. In addition, since z is AWGN, multiplying z by an 
orthogonal matrix does not change its statistics. It thus follows that the SER of S* and S are identical. 
To conclude the discussion of reduced constellations, an example of a constellation and its reduced 
constellation is provided. Consider S = {si,s 2 }, where Si = [vCU) \^5] T , and s 2 = [— Vol) — \/ol)} T . 
Using the definition of the reduced constellation, it is straightforward to see that the reduced constellation 
corresponding to S is the BPSK constellation set, and therefore the SER of S is identical to that of BPSK. 

IV. Symbol Error Rates of Multi-Dimensional Constellations 

Throughout the paper, the focus is to obtain a functional characterization of the SER of a multidi- 
mensional constellation, rather than to obtain bounds or closed-form expressions for the SER. Such a 
characterization can be used to uncover its convexity and complete monotonicity properties. 

A. Symbol Error Rates of Real Constellations Under AWGN 

To begin with, it is assumed the the transmitted symbol is a real vector, and the additive noise is white 
Gaussian. For constellations with reduced dimension N* = 1, the SER of the detector © under AWGN 
can be seen to be a positive linear combination of cm. functions of the form Q (J2pr]\ , rj > 0, which 
is cm.. The functional structure of the SER of constellations with N* > 2 is addressed in Theorem [TJ 
whose proof requires a result from the combinatorial geometry literature, which is stated next. 

Lemma 1. Let 7 be an N -dimensional polyhedron in Ji N with F facets. If CP contains the origin in its 



interior, then 7. admits the decomposition X := {Df q }f q , f = 1, ■ • • , F, q, — 1, . . . , Q ,, where T)f q 
are N -dimensional simplicial cones, and f can be viewed as an index of the facets of 7. 

Proof: An outline of the proof is provided in Appendix [A] ■ 
In other words, using the facets of an iV-dimensional polyhedron 7 C 7 N which contains the origin, 
it is possible to decompose 7 N into a collection of iV-dimensional simplicial cones. In what follows, the 
representation theorem is stated. 

Theorem 1. For a constellation S C 7 N , whose reduced constellation is S* and reduced dimension is 
N* > 2, the SER of the detector © under AWGN admits the representation 

P c (p) = fff cm (p), (11) 

where / cm (p) z * cm., and p > N*/2 — 1. In (ITTT) . the representing function of f cm (p) satisfies p,(u) = 
when u < d^jA, and jj,{u) > otherwise, where d m \ n is the minimum distance of the constellation. 

Proof: See Appendix [B] ■ 
To prove Theorem \T\ we work with a reduced constellation S* that is full rank. This is needed since 
Lemma [T]is used in Theorem [T] where 7 is a Voronoi region with dimension iV*. This highlights the need 
to work with the reduced constellation and dimension. The proof of Theorem Q] provided in Appendix |B] 
can be viewed as a generalization of the method adopted in [ETI to obtain an expression for the SER of 
arbitrary two-dimensional constellations under AWGN. 

As an example to corroborate Theorem [Q consider the square M -QAM constellation under AWGN. 
For this constellation, it is known that d m - m = a/2 and P e (p) — UiQ (y/rfp) — W2Q 2 (y/VP)^ w i m u i = 
4(VM- l)/\/M, tj = 3/(M - 1) and tu 2 = wf/4 (22||. In this case, P c (p) can be represented in the form 
(fTTI) . where p = (since A^* = 2), and the representing function of f cm (p) is given by 



*<(») = £ 

In 



'wiJ[0.5 < u < 1] (wi - w 2 )/[m > 1] 



(12) 



uy2u — 1 uy2u — 1 

which is zero when w < 0.5 and non-negative when u > 0.5, because U\ — w 2 > 0. Further, since p = 0, it 

is obvious that the square M-QAM constellation is an example of a constellation which has a cm. SER. 

According to Theorem \T\ the SER of every constellation can be written as P c (p) = pPfcmip), where 

p = N* /2 — 1 and / cm is cm.. However, it is possible that the SER of some constellations admit a 
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representation of the form ([TO . where the exponent of p is less than N*/2 — l. As a result, constellations 
for which this exponent is zero have cm. SERs. The following corollary of Theorem [1] establishes a 
necessary and sufficient condition for the SER of a constellation to be cm.. 

Corollary 1. The SER of § using the detector © under AWGN is cm. if its reduced dimension satisfies 
N* < 2. Conversely, let p,{u) = exp(pu)p~( N */ 2 ~^P e (p) dp, p{u) = J " p(u — v)v~ l ' 2 dv and r = 
\N*/2 — 2~|. If N* > 2 and N* is even (odd), the SER is cm. if and only if p(u) (p(u)) is r times 
differentiable, and d k p(u)/du k > (d k p(u)/du k > 0) for < k < r, and d r p(u) / du r (d r p(u) / du r ) is 
increasing and continuous. 

Proof: See Appendix O ■ 
In other words, Corollary [T] states that the complete monotonicity of the SER for constellations with 
N* < 2 does not depend on the geometry of the constellation. However, for constellations with higher 
reduced dimensions, the complete monotonicity of the SER depends on the differentiability of the rep- 
resenting function corresponding to f cm (p), which is a function of the constellation geometry and the 
a-priori probabilities. Although Corollary [T] is applicable to any constellation, it is not easy to obtain the 
equivalent set of conditions on the constellation geometry and prior probabilities under which the SER 
of a constellation with N* > 2 is cm., and this is posed as an open problem. 

Next, we provide instructive examples through which the complete monotonicity of the SER can be 
seen to depend on the constellation geometry for N* = 3. First, consider the constellation where the points 
are chosen as the vertices of a cube. In this case, the SER is given by P e (p) = 1 — (1 — Q (^y/2p)) 3 , which 
can be rewritten in the form ©, with p{u) = (3/7r)7[l < u < 2} + (ir - cos' 1 (a(u)))/2n 2 I[3 < u < 4] + 
(7t + cos _1 («(m)))/27t 2 /[m > 4], where a(u) = (3m 2 -12m + 8)/(m-2) 3 . Since p(u) > 0, from Bernstein's 
theorem, P e (p) is cm.. On the other hand, consider the 3-D square QAM constellation, whose points 
are given by all possible sign permutations of (±l/-\/6, ±l/\/6, ±l/\/6) and (±l/\/2, ±l/\/2, ±1/a/2). 
For this case, under the assumption of equal prior probabilities, numerical evaluation of the SER shows 
non-convexity (Fig. |3]). As a result, from ©, the SER is not cm.. Therefore, the cm. properties of the 
SER of a constellation with N* > 2 depends on the geometry and prior probabilities. 

With reference to existing literature, Corollary J] is a useful generalization of flU Theorem 2], which 
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does not address complete monotonicity, or the possibility of reduced dimension. Further, from Corollary 
[H it follows that the SER of any two-dimensional constellation under AWGN is convex. This particular 
consequence of Corollary [T] has been previously established in [[51 Theorem 1] using a different approach. 
In what follows, the behavior of second derivative of the SER for constellations with iV* > 2 is studied. 

Corollary 2. If the reduced dimension N* of a constellation S is greater than two, then the SER of the 
detector © under AWGN satisfies P " (p) > when p > p , where p := A(p + y^)/d^ in , p = N*/2 — 1, 
and d min is the minimum distance of the constellation. 

Proof: See Appendix [D] ■ 
A similar but weaker result has been obtained in [Q, where p = 4(iV + \/N) / c^ in . Corollary [2]provides 
a larger region where the SER has positive second derivative, which is an improvement over [5], since 
N*/2 - 1 < N* < N. Although Corollary [2] establishes a bound on the values of SNR for which the 
second derivative becomes non-negative, it does not forbid P e ' (p) > for p < p . Indeed, it is possible 
for some multi-dimensional constellation with iV* > 2 to posses a convex SER. An example of such a 
constellation is one with N = 3, wherein the points correspond to the vertices of regular convex poly tope 
(RCP) ifTTTl . for which convexity follows from the expression for the SER given by [fTTl Eqn. (2) - (7)]. 

B. Extension to Compound Gaussian Noise 

In what follows, Theorem[T]and its corollaries are generalized to the case of additive compound Gaussian 
noise. 

The system model considered is still as in ©, except that the additive noise is assumed to be z = y/W g, 
where W is a positive RV, which is independent of each component of g := [G\, . . . , Gn] t , and Gk ~ 
jV(0, 1/p) are iid, for k = 1, . . . , N. It should be noted that the elements of z are statistically dependent but 
uncorrelated in this case. Depending on the distribution of W, a number of noise distributions of interest 
arise from this formulation. For example, when W is an affine function of a Poisson RV, z follows a 
Middleton class-A distribution [fT3l . which is used to model multi-user interference; if W is a positively 
skewed alpha-stable RV with a characteristic function w( w ) = exp[— |u;| a (l — jsgn(aj) tan7ra/2)], where 
j = a/— T, < a < 1, and sgn(x) is the sign of x, then z follows the symmetric alpha-stable distribution 
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ffT4l . with characteristic function f z (uj) = exp[— |cu | 2ct ] . For different values of a, many impulsive noise 
distributions are obtained. For example, when a = 1/2, z is Cauchy distributed noise lfT4l . 

In this discussion, the minimum distance detector ©, which is still the maximum likelihood detector 
for dependent but uncorrected compound Gaussian noise, is assumed to be used at the receiver side. 
The reduced constellation corresponding to S is defined as in Definition Q] and is denoted by S*. It is not 
difficult to show that the SER of S and S* are identical under additive compound Gaussian noise. 

An extension of Theorem \T\ to the case of additive compound Gaussian noise is now developed. 
Conditioning on W = w, z is an i.i.d multivariate Gaussian, in which case Theorem Q] can be invoked to 
get P e (p\W = w) = fffcmip] w), where f cm (p; w) is cm. in p for each w. Averaging over the distribution 
of W, the equivalent of Theorem Q] for the case of compound Gaussian noise is obtained, since the 
expectation can be interpreted as a positive linear combination of f C m(p',w), which results in a cm. 
function, denoted by f cra {p)- 

Extensions of the corollaries of Theorem [T]for the case of compound Gaussian noise are also seen to 
be true, since they are obtained from the generalization of Theorem \T\ to this noise model, without any 
additional assumptions. 

C. Extension to Complex Constellations Under AWGN 

Theorem \T\ and its corollaries have been derived under the assumption that the transmitted symbol 
is chosen from a real constellation. A system model which is more relevant to in-phase/quadrature 
communication schemes and communication over fading channels is one where the transmitted symbol 
is a complex vector. Motivated by this, an extension of Theorem \T\ and its corollaries are obtained for 
the system model ©, where s is chosen from S = {si, . . . , s A/ } with E C N , i = 1, . . . , M, and the 
additive noise z ~ CA/"(0, (2/p)I). The receiver assumed for this discussion is the ML detector under 
circularly symmetric complex Gaussian noise, which is given by ©, where the ||-|| in © is interpreted as 
the 2-norm of a complex vector. Before proceeding to develop an extension of Theorem Q] to the complex 
case, a useful result pertaining to the SER of complex constellations is noted. 

Theorem 2. Let P Sji (p) denote the SER of © conditioned on the transmission ofsi e § := {s 1; . . . , s M }, 
where Sj e C N , j — 1,...,M, and z ~ £/V(0, (2/p)I). Further, let P e ,i(p) denote the SER of © 
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conditioned on the transmission of Sj G S := {§!,..., s M }, where Sj := [Re{Sj} T Im{s. ; } T ] T , j = 
1,...,M, and z ~ JV(0, (1//>)I). TTien, P M (p) = 

Proof: It is straightforward to show that the Voronoi regions corresponding to Sj and Sj are identical 
for i — 1, . . . , M. The Theorem then follows by observing that the SER is the Gaussian integral outside 
the Voronoi region. ■ 
Therefore, in order to generalize Theorem [T]to the case of a complex iV-dimensional transmitted symbol 
vector s G S, it is sufficient to consider the transmission of a real 2iV-dimensional vector s G S, where 
the i th element of S is given by [Re{s,} T Im{sj} T ] T , with G S for i = 1, . . . , M. Further, the SER S is 
equal to that of its reduced constellation S* described in Definition [fl and the reduced dimension N* of S 
is given by the rank of S, the i th column of which is given by [Re{sj} T Im{sj} T ] T , i = 1, . . . , M. Since 
§ is a real constellation, the expression for the SER is given by Theorem [H where N* is replaced by N*. 
From Theorem |2] this is also the SER of a complex iV-dimensional vector symbol set under complex 
AWGN. Using the same line of argument, it can be seen that the SER of a complex ^-dimensional 
constellation is cm. in the SNR if and only if its reduced dimension N* is one or two, as in Corollary 
[TJ Proceeding in a similar fashion, Corollary |2] can also be generalized to the complex case, by replacing 
N* with N* in the corollary to obtain the same conclusions. 

V. Applications 

In this section, applications of Theorem \T\ and Corollary [T]in the context of ordering of wireless system 
performance is presented. 

A. Applications in Stochastic Ordering 

Complete monotonicity of SER for complex constellations with a reduced dimension of one or two, as 
suggested by extension of Corollary [j] to complex constellations finds immediate application in comparing 
the average SER of such constellations over two different fading channels using the theory of stochastic 
ordering mentioned in Section III-BI To elucidate further, the following system model is considered: 

y = /is + v, (13) 
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where the effect of quasi-static fading is captured by the complex scalar RV h whose real and imaginary 
parts are independent of each other, s G § with s 6 C N , and z ~ CJ\f(Q, (l/p)I) is the circularly symmetric 
AWGN. For this system, the instantaneous channel gain is defined as X := \h\ 2 , and the instantaneous 
SNR is given by pX. Assuming that the receiver has full channel state information, the instantaneous 
SER can be shown to be a function of the instantaneous SNR only, and the average SER is obtained by 
taking the expectation of the instantaneous SER over the distribution of the instantaneous channel gain. 
In this application, goal is to compare the average SER of the system (fT3l) under two different fading 
channels with instantaneous channel gains X\ and X 2 . 

To begin with, let S be a complex constellation with a reduced dimension less than or equal to two, 
which has cm. SER according to the complex extension of Corollary \T\ Now, consider two fading scenarios 
with instantaneous channel gains Xi and X 2 , such that Xi < Lt X 2 . Then, according to ©, E [P c (pX 2 )} < 
E [P e (pXx)] , Vp > 0. As a result, complete monotonicity of SER can be exploited to compare two fading 
channels at all SNR, based on the average SER, even in cases where the expression for the average SER is 
not analytically tractable. As an illustrative example, consider the quadrature-PS K (QPSK) constellation, 
for which the reduced dimension can be seen to be equal to 2, and P e (px) = Q (y/2px) (assuming equal 
prior probabilities), which is cm. in x. Now, assume QPSK is used over two different Nakagami-m fading 
channels, the first one with LoS parameter mi and instantaneous channel gain X\, and the second one 
with LoS parameter m 2 and instantaneous channel gain X 2 , where m 2 > mi so that X\ <Lt X 2 . In this 
case, © implies that E [Q (y/2pX 2 )] < E [Q (y / 2pXf)] , Vp > 0, which provides a way of comparing 
the average SERs over two fading channels with different LoS parameters. 

The complex extension of Corollary [TJ suggests that the SER of a complex constellation with reduced 
dimension greater than or equal to four is not cm., and thus the LT ordering of instantaneous channel 
gains of two fading channels does not provide a conclusive comparison of the average SER of these 
channels. Motivated by this, a new stochastic order is introduced next, which can be used to compare the 
average SER of a multidimensional complex constellation over two different complex fading channels. 
The stochastic order <g p is formally defined below. 

Definition 2. Let X% and X 2 be two positive RVs, and let p > be fixed. Then Xi <g p X 2 if and only if 
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E [Xfexp(-pXx)] > E [X$exp(-pX 2 )] far all p > 0. 

In other words, for every p > 0, <g p is an integral stochastic order in the sense of @, with S P = 
{g(x) | g(x) = —x p exp(—px),p > 0}. A necessary and sufficient condition for X x < Sp X 2 can be proved 
to be as follows: 

Theorem 3. Let X\ and X 2 be two positive RVs, and p > 0. Then, X\ <g p X 2 if and only if 

E [Xff^X,)] > E [X 2 p / cm (X 2 )] , (14) 

where / cm (") z * c.m.. 

Proof: See Appendix El ■ 
Verifying the <g p order for any pair of random variables, when p E 3\fU{0} is relatively straightforward, 
and can be done by comparing the p th derivative of the real-valued Laplace transforms of the densities 
of the two RVs. Clearly, < Sp is the LT order, when p = 0. In this case, the envelope fading distributions 
for y/X~i and y/X^ such as Nakagami-m satisfy X\ <g X 2 , when mi < m 2 flU. Intriguingly however, 
for any p > 0, fading channels modelled using Nakagami-m or Rician distributions do not satisfy the 
< Sp order with respect to their corresponding line of sight parameters. For example, in the Nakagami-m 
fading scenario, X = \h\ 2 in (fT3l is Gamma distributed. In this case, 

E [** exp(-pX)] = m-{m + P )-^Y[m + p] 

T[m\ 

which increases with m for small p and decreases otherwise, for any fixed p > 0. Thus, if p > 0, 

Xi i % x 2 . 

Some implications of Theorem |3] to the ordering of average SERs of multidimensional constellations 
over fading channels are now considered. If X\ and X 2 are the instantaneous channel gains of two fading 
scenarios characterized by (fl3l . then according to Theorem [3l it is easy to show that 

X, < Sp X 2 E [P e (pXt)) > E [P c (pX 2 )) , Vp > , (16) 

where P e (•) is the instantaneous SER of a complex constellation S with reduced dimension N*, and 
p = N* /2 — 1. This is because, from the complex extension of Theorem Q] we have P c (p) = P p f C m(p), 
which implies P c (p) 6 {g(p)\g(p) = fffcmip)^ > }, and thus (fT6l) follows from Theorem |3] 
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It has been reported in the literature that, it is possible to find a fading distribution such that the SER 
of the AWGN channel is worse than that under the fading case at low SNR, when higher dimensional 
constellations are employed [5J. However, examples of such fading distributions have not been a subject 
of investigation. Using (fT6l ), it is now shown that the Nakagami-m fading case is an such an example. 
To begin with, consider the pure AWGN channel (i.e. the no fading scenario, where Pr[Xi = 1] = 1), 
for which E [Xf exp(— pXi)] = exp(— p). If X 2 denotes the instantaneous channel gain of a Nakagami- 
m fading channel, then for every < m < oo, we now argue that there exists a p\ > such that 
E[P c (pX 1 )] > E[P c (pX 2 )} for p < Pl , while E [P e {pXi)\ < E [P e (pX 2 )], for p > Pl . To this end, 
observe that E [Xf exp(— pXi)} is greater than E [X 2 exp(— pX 2 )} when p < p\, and vice-versa when 
P > Pi- Therefore, from (fl6l . the AWGN channel is worse than a Nakagami-m channel in terms of SER 
of constellations with N* > 2 at low SNR. 

Next, a relation between Xi < Sq X 2 and X\ < Sp X 2 is obtained, where p and q are non-negative. 

Theorem 4. Let X x and X 2 be two positive RVs. Then, for < q < p, X\ <g p X 2 X\ <g g X 2 . 

Proof: Since X x < Sp X 2 , from Theorem S we have E [Xff^Xx)} > E [X%f cm (X 2 )}, for every 
cm. function / cm (-). Choose f cm (x) := x^ k g cm (x), where < k < p, and g cm (x) as some cm. function. 
Clearly, f cra (x) as defined is cm., since x~ k is cm. for k > and a product of cm. functions is also 
cm.. As a result, according to Theorem |3] X% <g p X 2 implies E 
The theorem then follows by assuming q = p — k. ■ 
Theorem |4] in conjunction with (fTol implies that if X\ <g p X 2 , then X 2 is better than X\ in terms of 
average SERs of all constellations at all average SNR, with the reduced dimension of the constellation 
satisfying N* /2 — 1 < p. It is interesting to investigate the conditions on X\ and X 2 such that X\ <g p X 2 
for all p > 0. In that case, X 2 will be better than X x in terms of average SERs of any multi-dimensional 
constellation at all SNRs. However, this condition is not satisfied by any pair of random variables, as 
described in the following Theorem: 

Theorem 5. There are no two positive random variables which satisfy X\ <g p X 2 , for all p > 0. 

Proof: See Appendix IB ■ 
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As described in the paragraph above Theorem @] the AWGN channel is not the best channel in terms 
of SER of constellations with N* > 2 at all SNR. According to Theorem [5] there is no fading distribution 
which dominates every other fading distribution in the sense of <g p for all p > 0. Therefore, unlike cases 
where the SER metric is convex, where AWGN (no fading) outperforms any fading, this is not the case 
for N* > 2. Moreover, Theorem [5] suggests that there is no fading distribution which serves the role of 
"best" fading distribution, in terms of SERs of constellations of every dimension. 

VI. Conclusions 

In this paper, the SER of an arbitrary constellation with reduced dimension N* > 2 under AWGN 
is characterized as P c (p) = p p f C m(p)i where / cm (-) is a cm. function, and p > N*/2 — 1. This 
representation of the SER is shown to apply to cases when the noise follows a compound Gaussian 
distribution. The expression for the SER obtained herein is useful in establishing that the SER is a cm. 
function if the constellation has a reduced dimension of one or two. The complete monotonicity of SER 
for constellations with N* > 2 is shown to depend on the differentiability properties of the representing 
function corresponding to p~ p P c (p), which is a function of the constellation geometry and the prior 
probabilities. The exact relation between the constellation geometry and the complete monotonicity of 
the SER for constellations with iV* > 2 is left as an open problem. Complete monotonicity of the SER 
has applications in obtaining comparisons of averages of SERs over pairs of quasi-static fading channels, 
such as Nakagami-m, whose instantaneous SNRs are Laplace transform ordered. Such comparisons can 
be made even in cases where a closed-form expression for the average SER is not analytically tractable. In 
addition, a new stochastic ordering relation is introduced, which can be exploited to obtain comparisons 
of the average SER of an arbitrary multidimensional constellation over two different fading channels. 

Appendix A 
Proof of Lemma CD 

The proof of this lemma rests on the fact that it is possible to decompose Ji N into a set 3" consisting of 
polyhedral cones Ci, . . . , Cp, using the facets of an iV-dimensional polyhedron 7 G % N , which contains 
the origin in its interior [fT8l pp. 192]. Since every polyhedral cone admits a decomposition into N- 
dimensional simplicial cones |fT9l Lemma 1.40], it is possible to decompose each 6/ into a set of N- 
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dimensional simplicial cones {2)/,g } q , for / = 1, . . . , F. Consequently, Jl N admits a decomposition into 
iV-dimensional simplicial cones, given by • 

Appendix B 
Proof of Theorem CD 

Throughout this appendix, we work with the reduced constellation S*, with rank N*. Recall the AWGN 
system model ©. To obtain an expression for the symbol error rate averaged over all constellation points 
P e (p), we first evaluate P e>i (p) given by ®, and then use (fTOb . For the sake of simplicity, we assume 
that the Voronoi region of s* G S* is a polytope. The following proof can easily be extended to cases 
when %i is an unbounded polyhedron, by assuming an additional facet c x < 1, which turns %i into a 
polytope [18, pp. 75], and subsequently taking the limit of P c (p) so obtained as c — > 0. 

We begin with an outline of the proof. In general, evaluating © is not straightforward, since X$ is a 
polytope. However, the Gaussian integral in © can be simplified, if the region of integration is of the 
form Z — Z, where Z is an iV-dimensional simplicial cone, and Z is the intersection of a halfspace and Z. 
To this end, we show that the Voronoi region of s* has a dimension of N*, so that Lemma Q] can be used, 
and thus P e ,i (p) hi © can be rewritten as a sum of integrals over regions of the form Z — Z, where Z 
and Z are as defined above. Each of these integrals when expressed in hyperspherical coordinates yields 
a canonical structure, which can be algebraically manipulated to obtain (fTTT) in the Theorem. In order to 
show that N* /2 — 1 in (fTTT) is the smallest exponent for which the Theorem holds, an argument involving 
complete monotonicity of order a is provided towards the end of this appendix. 

In what follows, we present the details of the proof. Let %i be a non-redundant description of the 
Voronoi region of s* G §>*, for i = 1,...,M. With a slight abuse of notation, we are dropping the 
superscript * from UQ, to simplify the notation. First, we show that 3Q satisfies the conditions of Lemma 
[D To this end, for any set of affinely independent points in space (such as §>*), the dimension of the 
Voronoi region corresponding to each point is equal to the dimension of the afflne hull of the set (N*, 
when the set of points is S*) 11231 p. 232]. Therefore, %i is an N* -dimensional polytope in Jl N * . Also, 
since the origin of the coordinate system is shifted to s*, G Thus, satisfies the conditions of 
Lemma [Q As a result, using Lemma Q] we obtain a set Xj := {D f jq ,i} q j, which is a decomposition of 
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Jl N * into N* -dimensional simplicial cones (see Fig. 1). Clearly, every x ^ satisfies x G £>/,q,i — 
for some {g, /}. Let the number of facets of %i be F { . It now follows that 

Fi Qf,i 
/=1 9=1 

where 

J f^ = {£) N/2 J expf-^^jdn.-.d^. (18) 

In order to simplify the integral in ([TBI , we switch to the hyperspherical coordinate system [1241 . which 
is a generalization of the spherical coordinate system to higher dimensions. In this system, x G Jl N * is 
uniquely represented as [r, 1; . . . , ^» ], where r = ||x||, and 0^ is the angle between x and the k th 
edge of Df,q t i, k — 1, . . . , N* — 1. More precisely, let v fe define the unit vector in the direction of the 
k th edge of for jfc = 1, ... , N*. Then, fc = caT^x 1 ^,,/ ||x||), Jfc = 1, . . . , N* - 1 (See Fig. El for 

a two-dimensional example.). 

Next, we obtain the region of integration in ([TBI in hyperspherical coordinates. For any x G Df, q ,i — %i 
represented by [r, (j) x , . . . , JV ,_ 1 ], the parameter r must satisfy 

r f ,M) < r <°°> (19) 

where ff } qj(</)) is the distance of the point x G {x|aJjX = 6/,} from the origin. An expression for 
can t> e found by representing the hyperplane aj^x = 6/^ in hyperspherical coordinates, using 

the inverse hyperspherical transform relations 

fc-i 

x k = r cos 4> k Y[ sin 4> ki , k = 1, . . . , N* - 1 , 
fci=i 

2V*-1 

a; at* = r JJ sin fei , (20) 
fei=i 

and solving for r as a function of 4> := [0i, . . . , 0jv*-i]- Thus, we get 



Y, a h ,f,i cos (j) k [] s in fcl + a>N*-i,f,i H sin0 fei 



r /^(4>) = 3^=i —i ~ — > W 



k=l fci=l fcl=l 

In (l2~TT) . afcj^ is the fc th element of aT^ Also, for any x G — UQ, it is seen that 0^ must be 

at least radians (if x = av kq ,a > 0), and at most <Pk,f, q ,v which is the angle between v^, and 
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v kq , k — 1, . . . , N* — 1. In other words, 

< <t>u < cos" 1 (v^ g v fe , 9 ) =: 4 )/>9ii . (22) 

It is useful to note that 4>k,f, q ,i < 7r ' since it is the angle between any two edges of the simplicial cone 
D/ i9i i, which is at most %. 

Thus, (fT8l) can now be reformulated in hyperspherical coordinates, with the limits of integration given 
by (QS) and d22]) as 

4>N* -l,f,q,i 4>lJ,q,i OO 

J fM=(£) N ' /2 J ■■■ J J r^- 1 S (0)e-^ 2 drd0 1 ...d0 iVt _ 1 , (23) 

where s(0) := n^li" 2 sm^* - ^ -1 0^ arises from the Jacobian of the transformation. Substituting u = r 2 /2 
in (1231) . and changing the order of integration, we get 



P e ,, (p) = p N *' 2 J e-^jkM&u , (24) 

Q 

where p,i(u) is given by 



Ft Qf.i L 



>N*-l,f,q,i,l °l,f,q,i,l 



ftW:=^EEE / ... / sWvT^H [r f , q M? < «] • • • • (25) 

/=1 9=1 1=1 fl . J J 

¥-N*-l,f,q,i,l -l,f,q,i,l 

In (|25l) . L is the number of convex intervals of [<p\, . . . ,(j>N*-i] obtained after changing the order of 
integration, since the inverse function of F/ j(?i j(0) 2 is not unique. We now show that jli(u) > 0, which, 
together with Bernstein's Theorem implies that the integral in (124)) is equivalent to a cm. function of 
p. To this end, observe that the integrand in (12~5l is non-negative for u > 0, because s((f>) > for 
4>kj, q ,i G [0, 7r], A; = 1, . . . , N* — 1. Consequently, the result obtained after the iV* — 1 fold integration in 
(1231) is also non-negative. Thus, fii(u), which is a scaled version of a sum of non-negative integrals, is 
also non-negative. Therefore, through Bernstein's Theorem, we can assert that P e ,i (p) = p N *Wcm,i(p), 
where / cm ,i(p) is a c - m - function. 
Now, using (EH) in dTO]), we get 

OO 

P c (p)=p N */ 2 j e-^~p(u)du, (26) 
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where p(u) := Ei^i Pr [ s = s i]&( u )- Thus > (p) = p N * /2 fcm(p), where / cm (p) is cm. through 
Bernstein's Theorem, because /2(it) > as it is a positive linear combination of non-negative functions 

Pi(u),i = 1, . . . ,M. 

Next, we strengthen the representation (|26l by showing that p a f C m(p) is c - m - for a = 1. To this end, 
recall from Section Hi- A I that the necessary and sufficient condition for a cm. function to be cm. of order 
1 is that its representing function be nonnegative and increasing. This is indeed the case for jl(u). Thus, 
we have just showed that f cm (p) is cm. of order a — 1. Denoting f cm (p) := pf C m(p), which we have just 
showed to be cm., we get a stronger representation for the SER using (T2~6l) as follows 

iV* i 

P c (p)=p-- 1 f cm (p) i (27) 

where f cm (p) is cm.. 

Next, the support of the representing function of / cm (p) is investigated. Let p(u) be the representing 
function of f C m(p)- Accordingly, by applying integration by parts on (|26l ). it is seen that 

M E, Qf,i L 

m«) = EEEE / ••• / [^w 8 <«] + 

i=l /=1 g=l J=l fl J fl J 

S.N*-l,f,q,i,l 2.1, f,q,i,l 

u^- l I [f Uq Mf = u])d<P 1 ... d^,^ , (28) 

which is zero if u < min inf Tf q ^(4>) 2 , and non-negative otherwise. Recalling the expression for ff, q ,i(<p) 
from (|2TI) . it is immediately seen that min inf Tf qi ((f>) 2 = min/j b 2 fi . Observing that min/j&^ = d^jA, 

q,f,i <f> ■>' J > 

we conclude that the support of [i is contained in [^ in /4, oo). 
This concludes the proof of the Theorem. 

Appendix C 
Proof of Corollary CD 

The direct part of the corollary for constellations with N* = 1 is immediate, since P c (p) is a positive 
linear combination of functions of the form Q ( v /2pr/) ,rj > 0, which is known to be cm. (see e.g. JH). 
Also, the complete monotonicity of the SER for the case of iV* = 2 is straightforward from Theorem [TJ 
To see the converse, we use the necessary and sufficient condition for a function f(p) to be cm. of order 
a, as described in Section Hl-AI For the case when iV* is even, we assume f(p) = f cm (p), where / cm is 
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as defined in Theorem [U and a = N*/2 — 1. For the case of odd N*, we assume f(p) = p 1 ^ 2 f C m(p), 
and a = \N*/2 — 1]. The converse thus easily follows. 

Appendix D 
Proof of Theorem [2] 

Recall from Theorem Q] that for a constellation with reduced dimension iV*, the SER can be written as 

oo 

Pe (p)=P P J ex V (-pu)p(u)du , (29) 



where p = N*/2 — 1, with p{u) being zero when u 6 [0, c^ in /4), and non-negative otherwise. We now 
obtain sufficient conditions for P e (p) > 0. Differentiating (|29l twice with respect to p, we obtain 

oo 

P e " (P) = P p - 2 J e"">(«) (r* - P -^-) (u - P -^j du , (30) 



where the differentiation under the integral sign in (|29l is permitted, because the limits of integration 
are independent of p. A sufficient condition for (|30l) to be non-negative is that the integrand in (l30l is 
non-negative. Accordingly, observe that the integrand in (|30l) is non-negative when u > (p + ^/p) / p. 
Thus, if the support of p is a subset of \{p + y/p)/p, oo), it follows that P " (p) > 0. In other words, if 
p > 4(p + v /p)/c?min' we have P e (p) > 0, which proves the result. 

Appendix E 
Proof of Theorem [3] 

For any two non-negative RVs X x and X 2 , we have 

< % X 2 & E [X{ exp(-pX 1 )] > E [X p 2 exp(-pX 2 )] , Vp > 0. (31) 
We now establish the following: 

X x < Sp X 2 ^E [X p J cm (pX 1 )} > E [X p 1 f cm (pX 1 )} Vp > 0, (32) 
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where / cm (p) := / °° exp(— pu)p(u)du is a cm. function, for some p(u) > 0. Assume u > 0, then 
Xi < 9p X 2 E [Xf exp(-pXiu)] > E [X% exp(-pX 2 u)} Vp > 0. Next, observe that 



E [Xff^pX^] = E 



oo oo 

J Xfexp(-upX 1 )^(u)du =y > E[X 1 7 cm (pX 1 )]p( M )d M (33) 



oo 

> y E [Xf/ cm (pX 2 )] p( M )d M = E [Xf/ cm (pX 2 )] , (34) 



o 



Vp > 0, provided the expectations exist. This proves the direct part of the theorem. To see the converse, 
let /cm(p) = exp(— px), which is cm. in p for each x > 0. 

Appendix F 
Proof of Theorem [5] 

Let Xi, X 2 be non-negative RVs with PDFs f Xl (x) and fx 2 (%) respectively. We now show that 
X\ <g p X 2 cannot hold for all p > 0, by showing that X\ <g p X 2 does not hold for every p in a subset 
or Jl U {0}, i.e., for p G N U {0}. In order to satisfy Theorem [3] for every p G N U {0}, we require for 
every p 

(-l)P— [E [exp(-pXa)] - E [exp(-pX 2 )]] > Vp, (35) 

where we have used the identity J °° s p exp(— ps)p(s)ds = (—l) p (d p /dp p ) / °° exp(— ps)p(s)ds. Recalling 
the definition of a cm. function from ([T]), we gather that E [exp(— pXi)] — E [exp(— pX 2 )} in (1351) must 
be a cm. function. Consequently, Bernstein's Theorem mandates that f Xl (x) — fx 2 ( x ) > Vr. However, 
this condition is never satisfied by any pair of random variables, since both the density functions must 
individually integrate to unity, which cannot be the case if f Xl (x) > fx 2 (x) Vx. Thus, the Theorem 
follows. 
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Fig. 1. Voronoi region for signal point s* £ §*, where the reduced dimension of the constellation S is N* = 2. The origin of the coordinate 
axes is shifted to s*. Using Lemma 1, 3i N is decomposed into a collection of 2-dimensional simplicial cones Xj := . . . , T>6,i,i} 

using the facets of 3Q. 




Fig. 2. The 2-dimensional simplicial cone Di^.i obtained using Lemma 1 with Xj, represented in hyperspherical coordinates. 
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Fig. 3. Second derivative of P c (p) for the 3-dimensional square QAM constellation, 
constellation with TV* = 3 whose SER is not cm.. 



Since the SER is not convex, this is an example of a 



